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^ : Abstract 

We propose the minimal supersymmetric sterile neutrino model (MSSNM) where the sterile 
neutrino masses are about 1 eV, while the active neutrino masses and the mixings among the active 

o 

lO ■ and sterile neutrinos are generated during late time phase transition. All the current experimental 

o 

neutrino data include the LSND can be explained simultaneously, and the constraints on the sterile 
neutrinos from the big bang nucleosynthesis and large scale structure can be evaded. To realize 



^ ' the MSSNM naturally, we consider the supersymmetric intermediate-scale U{iy model, the low 

energy C/(l)' model with a secluded C/(l) '-breaking sector, and the DFSZ and KSVZ axion models. 
• i-^ ■ 

^ I In these models, the /x problem can be solved elegantly, and the 1 eV sterile neutrino masses can 

be generated via high-dimensional operators. For the low energy model with a secluded 

C/(l) '-breaking sector, we also present a scenario in which the masses and mixings for the active 
and sterile neutrinos are all generated during late time phase transition. 

PACS: 12.60.Jv; 12.60.Cn; 14.80.Mz 
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I. INTRODUCTION 

There has been great progress in neutrino physics during last several years Solar 
neutrino and atmospheric neutrino [3| experiments together with reactor neutrino 
experiments have established the oscillation solutions to the solar and atmospheric neutrino 
anomalies, which are consistent with three light active neutrino scheme. However, the LSND 
experiment found evidence for the oscillations —>■ Ue and —>■ z/g with an oscillation 
probability of around 3 x 10"'^ [6] and a Am^ ^ leV^. The statistical evidence for the anti- 
neutrino oscillations is much stronger than that for the neutrino case, with some analyses 
finding a 5cr effect f^. Although the other experiments eliminated a large fraction of the 
parameter space allowed by the LSND, they do not exclude the LSND result ^. Hopefully, 
the Mini-BOONE experiment at Fermilab will settle this issue down in the near future 9|. 

If the LSND experiment is confirmed, to explain its data, we need introduce one or two 
sterile neutrinos, which are Standard Model (SM) singlets and can mix with the active neu- 
trinos. Also, the masses for the sterile neutrinos should be in the eV range. Because of 
various possible mass hierarchies for the active and sterile neutrinos, there are three propos- 
als: the 2+2 model (ld|, 3+1 model 

seems less disfavored than the 2+2 model due to the null results of other oscillation experi- 



12| is apparently in better agreement with all the experimental 



ments. And the 3+2 model 
data than the others. 

There are two strong constraints on the sterile neutrino models which can explain the 
LSND result. The first constraint is that the big bang nucleosynthesis (BBN) allows about 
three effective light neutrinos {N'^^^) in the equilibrium when the Universe temperature is 
around 1 MeV [13]. However, for above three proposals, the rapid active neutrinos-sterile 
neutrino(s) oscillations give the N^ff = 4 for the 3+1 and 2+2 models, and Nf^ = 5 for 
the 3+2 model. The second constraint is the bound on the sum of all the neutrino masses 
in the equilibrium at the epoch of structure formation which corresponds to a temperature 
around an eV from the large scale structure surveys and WMAP [lj,13- Suppose that the 
sterile neutrinos are in the equilibrium during the BBN epoch, the upper bound on the sum 
of all the neutrino masses is about 1.38 eV for the 3+1 and 2+2 models, and about 2.12 
eV for the 3+2 model. These constraints are also very severe because for example the 3+2 
model is close to be ruled out if we take the face value. 
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To avoid the BBN and large scale structure constraints and explain the LSND experiment, 
Chacko, Hall, Oliver and Perelstein proposed a class of models where the masses and mixings 
of the active and sterile neutrinos are generated during late time phase transition Based 
on the next to the minimal supersymmetric Standard Model (NMSSM), they introduced two 
SM singlet fields with vacuum expectation values (VEVs) in the 100 keV range so that at 
the BBN epoch the active as well as the sterile neutrinos are massless. Thus, there is 
no oscillation among them which can bring the sterile neutrinos into equilibrium. Since the 
sterile neutrinos decouple from Hubble expansion at very high temperatures, their abundance 
at the BBN epoch is suppressed leading to concordance with the BBN constraints. Moreover, 
the constraints on the sum of the neutrino masses from large scale structure surveys and 
WMAP are also easily avoided. The breaking of the global symmetries gives rise to a few 
Goldstone bosons, which can couple to both active and sterile neutrinos. These couplings 
are strong enough for the sterile neutrinos to disappear after they become non-relativistic, 
for example, by decaying into an active neutrino and a Goldstone boson. Thus, the relic 
abundance of the sterile neutrinos is low, and they do not significantly contribute to dark 
matter. 

Furthermore, using the idea of the late time phase transition, Mohapatra and Nasri 
considered the sterile neutrinos in the mirror matter models R. If the sterile neutrinos are 
the mirror neutrinos, they only need to generate the mixings between the active and sterile 
neutrinos (and not their masses) via the late time phase transition to avoid the constraints 
from the BBN, the large scale structure surveys and WMAP. An advantage of this model 
is that the contribution of the sterile neutrinos to the energy density of the universe at 
the BBN epoch is given by a free parameter unlike the model in Ref. 16|. Cosmological 
consequences of these two kinds of the models have also been discussed in Refs. 

00. 

In this paper, we propose the minimal supersymmetric sterile neutrino model (MSSNM) 
with late time phase transition. The active neutrino masses and the mixings among the 
active and sterile neutrinos are generated during late time phase transition, while the 1 eV 
masses for the sterile neutrinos are introduced directly in the Lagrangian. We can also forbid 
the dangerous operators by introducing Z3 x Z2 global symmetry. In the MSSNM, the current 
neutrino data from all the experiments include the LSND can be explained simultaneously, 
and one can automatically evade the constraints on sterile neutrinos from the BBN, large 
scale structure surveys and WMAP. However, there are two interesting questions in the 
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MSSNM: (1) how to produce the 1 eV masses for the sterile neutrinos? (2) how to solve the n 
problem because the supersjTiimetry breaking is mediated by gauge interactions? To realize 
the MSSNM naturally, we consider the supersymmetric intermediate-scale U{1)' model, 
the supersymmetric low energy f/(l)' model with a secluded f/(l)'-breaking sector, and the 
supersymmetric Dine-Fischler-Srednicki-Zhitnitskii (DFSZ) and Kim-Shifman-Vainshtein- 
Zakharov (KSVZ) axion models Q|. In these models, the 1 eV masses for the sterile 
neutrinos can be obtained via the high- dimensional operators by integrating out the heavy 
fields, and the dimension-5 operators for the active neutrino masses and the mixings among 
the active and sterile neutrinos can also be generated by integrating out the heavy fields. 
Also, the fi problem can be solved elegantly. Furthermore, for the low energy model 
with a secluded U (l)'-breaking sector, we briefiy present a scenario where the sterile neutrino 
masses are also generated during late time phase transition. 

This paper is organized as follows: in Section II, we propose the MSSNM. We consider 
the supersymmetric intermediate-scale f/(l)' model, the supersymmetric low energy f/(l)' 
models with a secluded f/(l)'-breaking sector, and the supersymmetric DFSZ and KSVZ 
axion models in Sections III, IV and V, respectively. Our discussions and conclusions are 
given in Section VI. 

II. MINIMAL SUPERSYMMETRIC STERILE NEUTRINO MODEL WITH 
LATE TIME PHASE TRANSITION 

We first specify our conventions. For the supersymmetric Standard Model, the SM 
fermions and Higgs fields are superfields belonging to chiral multiplets. The left-handed 
quark doublets, the right-handed up-type quarks, the right-handed down-type quarks, the 
left-handed lepton doublets, the right-handed neutrinos, the right-handed leptons, and one 
pair of Higgs doublets are denoted as Qi, u^, d^, Li, rii, e^, Hu and Ha, respectively. 

To construct the MSSNM with late time phase transition, we introduce one SM singlet 
field 0. The relevant superpotential is 

(J) K n 

= XijLiTijHu— + -0^ + mnijTiinj , (1) 

where Ajj and k are the Yukawa coupling constants, and m„jj are the masses for the ster- 
ile neutrinos that are around 1 eV and can be generated via high-dimensional operators 



4 



after extra gauge or global symmetry breaking in the following model buildings. The non- 
renormalizable term (the first term) can be obtained at scale M by integrating out the heavy 
fields where M is around 10^"^ GeV. 

After the electroweak symmetry breaking, the renormalizable effective Lagrangian for the 
neutrino sector is 

- Ly = {gijUiUjcj) + mnijUiUj + H.C.) + V(0) , (2) 

where gij = Xij < > /M, z/j is the left-handed neutrino, and the scalar potential is 
V(0) = — + Here, we have assumed that the supersymmetry breaking effects 

produce the negative soft mass-squared for 0. Thus, there is one global U{1) symmetry 
under which z/j and have opposite charges with the same magnitude while nj is neutral. 
When (j) acquires VEV, this global U{1) symmetry is broken. And then, the active neutrinos 
obtain masses, and there is one pseudo-Goldstone boson which has diagonal couplings to 
the neutrinos in the mass basis. 

Similar to the discussions in Ref. we can avoid the constraints from the BBN, the 
large scale structure surveys and WMAP. First, because the total energy density in radiation 
at the time of BBN does not differ significantly from the SM prediction we require that 
the "hidden sector" fields Ui and not be in thermal equilibrium with the "observable 
sector" fields (z/j and 7, etc) before and during the BBN. To be concrete, we require that the 
two sectors decouple at a certain temperature Tq > 1 GeV, and do not recouple until the 
temperature of the observable sector drops below Tw ~ 1 MeV, the temperature at which 
the weak interactions decouple. Thus, we have 

gij < 10-5 , gijK < lO-^V-i , (3) 

where r is the ratio of temperature of the hidden sector to that of the observable sector at 
the time of BBN. The energy density in the hidden sector is suppressed by a factor of 
compared to that of observable sector from the naive estimation, so, r ^ 0.3 is enough for 
one to avoid the BBN constraints. Moreover, because at least one active neutrino has mass 
around 0.05 eV, we obtain that / ^ 10 keV. And to avoid producing sterile neutrinos by 
oscillations prior to the decouplings of weak interactions, we obtain that / ^ r MeV, and 
then g-ij ^ r^^ 10^'' and k, ^ 10^'^. In short, from the BBN constraints, we have 

10 keV ^ / ^ rMeV , 10"^ ^ gij ^ 10"^ , gtjK ^ IQ-^^r-^ . (4) 
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Second, the constraints on the sum of neutrino masses from the large scale structure 
.u„ey. and WMAP Q a. autcatically .vaded in above .ode., and do no. lead to e... 
limits on /. The above lower bound on gij implies that the reactions uu ^ nn, 00 become 
unfrozen before the sterile neutrinos become non-relativistic. These reactions thermalize the 
hidden sector fields with the active neutrinos. The density of the thermal sterile neutrino 
with mass at temperatures T < is suppressed by a Boltzmann factor e""^"/"^, and the 
excess sterile neutrinos disappear either via a decay process n z/0, or via an annihilation 
process nn vv. Thus, the massive sterile neutrinos will not give a significant contribution 
to dark matter. And only the sum of the masses of the active neutrinos and the Goldstone 
boson has to satisfy the constraints in Ref. jisj ]. 

In our model, the supersymmetry breaking must be mediated via the gauge interactions 
jioj l. And only feels the supersymmetry breaking via its coupling to Li and so that its 
supersymmetry breaking soft mass is around 100 keV. However, the /i problem becomes a 
severe problem because of the gauge mediated supersymmetry breaking. 

In addition, in the MSSNM, we must highly suppress some other renormalizable operators 
in the superpotential which are allowed by the gauge symmetry, for example, 0^, 0^nj and 
if„Lj0, etc. To achieve this, we introduce a global Z3 x discrete symmetry. Under the 
Z3 symmetry, the particles in the MSSNM transform as 

(Q„ n„ (Q„ n„ , «, Hi) e-^^./a (^c^ ^ 

(0, d% H^) e^2-/3(0, d'r, L„ H^) . (5) 

And under the Z2 symmetry, the particles in the MSSNM transform as 

(Qi, <, rf., L,, rii, ef) — > - {Qi, m-, d% n^, e-) 

(0, Hu, Hd) (0, Hu, Hd) . (6) 

Note that the jiHuHd term in the superpotential is allowed by the x Z2 discrete symmetry. 

To solve the problem and generate the 1 eV sterile neutrino masses naturally via the 
high-dimensional operators, we shall consider the extra gauge symmetry or global symmetry. 
In the following Sections, we will show that the MSSNM can be realized elegantly in the 
intermediate-scale U{iy model, the low energy f/(l)' model with a secluded f/(l)'-breaking 
sector, and the axion models. In these models, we do not need to introduce above Z^ x Z2 
discrete symmetry. 
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TABLE I: The {/(I)' charges of the relevant particles in the intermediate-scale 17(1)' model. Here 
Qs and Ql are the f7(l)' charges for S and Li, respectively. 



Field 




S' 


</> 


Hu 




X 


X 


X' 


x' 


Charge 




-3Qs 




-Ql + 17Q5/6 


Ql - 29Q5/6 


-nQs/Q 


IIQ5/6 


-2Qs/3 


2Qs/3 



III. INTERMEDIATE-SCALE [/(I)' MODEL 

We first consider the intermediate-scale U{iy model where the /i problem can be solved 
simultaneously. To break the f/(l)' gauge symmetry, we introduce two SM singlet Higgs 
fields S and S' with f/(l)' charges Qs and —3Qs, respectively. After the supersymmetry is 
broken, the Higgs potential for 5* and S' is 

V{S, S') = ml\S\' + m%\Sf + i^|,g| {\S\' - 3|5f , (7) 

where m| and m|, are the supersymmetry breaking soft masses for S and S', respectively. 
To break the intermediate-scale f/(l)' gauge symmetry, we assume that the sum of the 
supersymmetry breaking soft masses for S and S' is negative, i.e., m| + m|; < 0. Then, 
there is a runaway direction along the D-flat direction -\/3K5')| = |(5")|. However, the 
potential can be stabilized by the loop corrections or higher-dimensional operators. Thus, 
the S and S' fields can acquire the intermediate-scale VEVs. For example, suppose that 
there is a high- dimensional operator in the superpotential 

where Mpi is the Planck scale, we obtain that the S and 5" fields can have the VEVs around 
10^° GeV. So, the U{1)' gauge sjTumetry is broken at intermediate scale elegantly. 

To realize the MSSNM naturally, we introduce four SM singlet fields X, X, X', and x' . 
The U{iy charges for the relevant particles are given in Table HI 

The relevant superpotential is 

W = y,L,H^X + y'^n,<f)X + yxSXX + yx'X'S<f) + yj^>X'<j)^ + Mx'X'X' 



where Ui, y'j, yx, Vx', Vx', Kj and h are Yukawa couplings, and Mx' is the vector- hke mass 
for X' and x' . We assume that Mx' is about 10^^ GeV, which can be generated from the 
non-renormahzable operators after the gauge symmetry in the Grand Unified Theory (GUT) 
is broken because Mx' ~ MQ^rp/ Mpi where Mqut is the GUT scale. 

After the f/(l)' gauge symmetry breaking, we obtain the superpotential 

W = XijLiTljHuYY- + n^^ + ^nijUiTlj + jJ^HuHd , (lO) 
IVlx o 



where 



(S) 

\j = ViVj , Mx = yx{S) , K = Syxfyj^'—— , (11) 

Mx' 



{sr ,. AS)' 



\2 

-3 



To produce the suitable sterile neutrino masses rrinij, we need Xij ~ 10 . Such a value 
for Xij could be generated if the corresponding operators for sterile neutrino masses in Eq. 



were themselves due to the high-dimensional operators involving additional 



fields with 



2^. 



VEVs close to Mpi, e.g., associated with an anomalous f/(l)' gauge symmetry 

We do not consider the U{1)' anomaly cancellation in this paper because the anomaly 
free U{1)' models can be constructed easily by introducing SM vector-like fields if one follows 



the procedures in Refs. |2lL l22|. 



IV. LOW ENERGY C/(l)' MODELS WITH A SECLUDED f/(l)'-BREAKING SEC- 
TOR 

In the low energy supersymmetric U{iy models, the fi problem can be solved elegantly 
with an effective fi parameter generated by the VEV of the SM singlet field 5* which breaks 
the U{iy symmetry. And the Minimal Supersymmetric Standard Model (MSSM) upper 
bound of Mz on the tree-level mass of the corresponding lightest MSSM Higgs scalar is 
relaxed because of the Yukawa term hSHdHu in the superpotential and the U{iy D- 
term [24 1. More generally, for specific f/(l)' charge assignments for the ordinary and exotic 
fields one can simultaneously ensure the absence of anomalies; that all fields of the low energy 
effective theory are chiral, avoiding a generalized /i problem; and the absence of dimension-4 



proton decay operators 



2l|. 



There are stringent limits from direct searches at the Tevatron 23 and from indirect 
p.eo.i„n tes. at .he Z-po.e. a. LEP 2. a.d .o. weaU „eut.a. cu„e„t experiments Q. 

The constraints depend on the particular Z' couplings, but in typical models one requires 
Mz' > (500 — 800) GeV and the Z — Z' mixing angle az-z' to be smaller than a few x 10^'^. 
To explain the Z — Z' mass hierarchy, Erler, Langacker and Li proposed a supersymmetric 
model with a string-motivated secluded f/(l)'-breaking sector, where the squark and slepton 
spectra can mimic those of the MSSM, the electroweak symmetry breaking is driven by 
relatively large A terms, and a large Z' mass can be generated by the VEVs of additional 
SM singlet fields that are charged under the U{1) The phenomenological consequences 
of the low energy f/(l)' models, especially the models with a secluded f/(l)'-breaking sector, 



have been studied extensively 
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First, let us briefly review the supersymmetric f/(l)' model with a secluded f/(l)'-breaking 
sector j27|. There are one pair of Higgs doublets Hu and Hd, and four SM singlets, S, Si, 
S2, and 5*3. The f/(l)' charges for the Higgs fields satisfy 

Qs = -Qs, = -Qs, = IQss , Qh, +Qh^+Qs = 0. (13) 

The superpotential for the Higgs fields is 

Wh = hSH^Hu + , (14) 

where the Yukawa couplings h and A are respectively associated with the effective // term 
and with the runaway direction. The corresponding F-term scalar potential is 

Vj. = {\Hd\W + \s\W + \s\W) 

+X' {\Si\'\S2\' + \S2\'\Ss\' + l^an^ip) , (15) 
And the D-term scalar potential is 



(16) 



where G"^ = + g"^; gi,g2, and gz' are the coupling constants for U{1), SU{2)i and U{1)' , 
respectively; and is the U{1)' charge of the field 
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In addition, we introduce the supersymmetry breaking soft terms 

3 

Kf/t = mjjjHd\^ + mjjjH^\' + ml\S\^ + J2m%\S,\^-{AhhSHdH^ + AxXSiS2S3 

1=1 

+mls,SS, + mls,SS2 + li.C.) . (17) 



There is an almost F and D flat direction involving Si, with the flatness lifted by a small 
Yukawa coupling A. For a sufficiently small value of A, the Z' mass can be arbitrarily large. 
For example, if /i ~ lOA, one can generate the Z — Z' mass hierarchy where the Z' mass 
{Mz') is at the order of 1 TeV 

We shall consider the model where the sterile neutrinos are charged under the f/(l)' 
gauge symmetry. Because the sterile neutrinos must be decouple with the observable sector 
fields at Universe temperature about 1 GeV or at least before the QCD phase transition 
with chiral symmetry breaking, the U{1)' gauge interaction must be decouple at the same 
temperature. For the f/(l)' interaction mediated by massless gauge boson, the interaction 
rate F ~ n(7\v\ ~ o'z'T where az' = gz'/{4:'n'). And during radiation-dominated epoch, 
H ~ T^/Mpi. So, T/H ~ al,Mpi/T, and for T > 10^^ GeV such reaction is effectively 
decouple. For the U{1)' interaction mediated by massive gauge boson, the interaction rate 
F ~ G|,T5 where Gz' = az'/M^,. And then T/H ^ G%,MpiT\ Thus, if we require that 
the U{1)' gauge interaction decouple at temperature about 1 GeV, 

Tdecoupie ~ [j^^^) MeV > 1 GeV , (18) 

we obtain that the f/(l)' gauge boson mass Mz' is larger than about 18 TeV. And if Tdecoupie > 
300 MeV, we find that Mz' > 7.2 TeV. These roughly estimations agree with the relevant 
detail calculations in Ref. j29|. 

In the U{iy model with a secluded [/(I) '-breaking sector, by running our code, we have 
shown that we can indeed generate the 10 TeV scale mass for the f/(l)' gauge boson by 
choosing smaller A and suitable supersymmetry breaking soft parameters. Let us give an 



exmple. We choose the standard GUT value gz' = y5/3gi (It is y5/3(7i that unifies with 
g2 and gs in the simple GUT models.). With the input parameters in Tab el ITTj we obtain 
the VEVs for the Higgs fields after minimizing the Higgs potential numerically: {H^l) = 123 
GeV, {H^) = 123 GeV, (S) = 128 GeV, {Si) = 10846 GeV, (^2) = 10846 GeV, (^3) = 
10846 GeV. And then, we get that the f/(l)' gauge boson mass is 17334 GeV. 
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TABLE II: The input parameters and Higgs VEVs in the 11(1)' model with a secluded 
breaking sector. For the mass parameters and the mass-squared parameters, the units are GeV 
and GeV^, respectively. 



h 


A 


Ah 


Ax 


2 

rriTT 




0.9 


0.01 


219 


219 


-15.52 


-15.52 
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'"is. 


-21.92 


19.32 


19.32 


-10.92 


2.42 


2.42 






Qs 


Qs, 


Qs2 


Qs, 


0.5 


0.5 


-1 


1 


1 


-2 




im) 


{S) 


{Si) 


{S2) 


(53) 


123 


123 


128 


10846 


10846 


10846 



In addition, if (p is charged under the U {!)' gauge symmetry, its supersymmetry breaking 
soft mass will be much larger than 100 keV without fine-tuning due to the low energy C/(l)' 
gauge interaction. Thus, unlike the intermediate-scale 11(1)' model, we assume that (p is 
neutral under the 11(1)'. 

To realize the MSSNM, we define a global symmetry in our model 

m^m, e-''^/'Si , $ e'^'^/^^ , (19) 

where $ denotes all the other fields. One can easily check that the above superpotential and 
supersymmetry breaking soft terms satisfy this Z3 global symmetry. 
With the following 17(1)' charges for the relevant fields, 

Qn, = - IQs, , Qh^ + Ql, = IQs, , (20) 

we can have the superpotential for neutrino sector 

= K,L,n,H^^ + i*^' + i (^i*^' + ^^-^2) n^n, , (21) 

where the non-renormalizable operators can be obtained by integrating out the heavy fields 
at scales M and M'. Here, M M' 10^~^ GeV. Thus, after the U{1)' gauge symmetry is 
broken, we obtain 

= ^{^us^)'+m^^)') ■ (22) 
11 



Furthermore, we can construct the model where the sterile neutrino masses are also 
generated during late time phase transition. For example, with the following U{iy charges 
for the relevant fields, 

Qni = - 'jQsi , Qhu + Ql, = 2^Si , (23) 
we can have the superpotential for neutrino sector 

The discussions for this model are quite similar to those in Section II, so, we will not give 
them here. 



V. DFSZ AND KSVZ AXION MODELS 



As we know, the strong CP problem is solved elegantly by the Peccei-Quinn (PQ) mecha- 
nism js^, in which a global axial symmetry U{1)pq is introduced and broken spontaneously 
at some high energy scale. The original Weinberg-Wilczek axion 31] is excluded by experi- 
ment, in particular by the non-observation of the rare decayi^-^vr + a[3l. And there are 
two viable "invisible" axion models in which the experimental bounds can be evaded: (1) the 
DFSZ axion model, in which a SM singlet and one pair of Higgs doublets are introduced, 
and the SM fermions and Higgs fields are charged under U{1)pq symmetry (2) the 
KSVZ axion model, which introduces a SM singlet and a pair of extra vector-like quarks 
that carry t/(l)pQ charges while the SM fermions and Higgs fields are neutral under U{1)pq 
symmetry 19|. In addition, from laboratory, astrophysics, and cosmology constraints, the 
U{l)pQ symmetry breaking scale fa is limited to the range 10^° GeV < fa < 10^^ GeV [3^ . 

The quantum gravitational effects, associated with black holes, worm holes, etc., are be- 
lieved to violate all the global symmetries, while they respect all the gauge symmetries . 
These effects may destabilize the axion solutions to the strong CP problem due to the 
violation of the global Peccei-Quinn symmetry. However, after a gauge symmetry is spon- 
taneously broken, there may exist a remnant discrete gauge symmetry which will not be 
violated by quantum gravity (3J]. Thus, we can avoid the destabilization problem asso- 
ciated with quantum gravity by introducing an additional approximate global symmetry 
arising from the broken gauge symmetry. 
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In string model buildings, there generically exists at least one anomalous U{1)a gauge 
symmetry with its anomalies cancelled by the Green-Schwarz mechanism js^. The anoma- 
lous U{1)a gauge symmetry is broken near the string scale when some scalar fields, which are 
charged under U{1)a, obtain VEVs and cancel the Fayet-Iliopoulos term oiUil)^- Then 
the D-fiatness for U{1)a is preserved and the supersymmetry is unbroken jsGj. Usually, 
there is an unbroken discrete Zn subgroup of the U{1)a gauge symmetry, which is protected 
against quantum gravitational violation. We shall consider this Zn discrete symmetry as 
an additional global symmetry to forbid the dangerous non-renormalizable operators which 
can destabilize the axion solutions to the strong CP problem 

For the gauge symmetry JJiGi x U{1)a, the Green-Schwarz anomaly cancellation condi- 
tions from an effective theory point of view are Isol . I40I 1 



^gravity 



Sgs , (25) 



ki 12 

where the Ai are anomaly coefficients associated with Gf x U{1)a, ki is the level of the 
corresponding Kac-Moody algebra, and 6gs is a constant which is not specified by low- 
energy theory alone. For a non-Abelian group, ki is a positive integer, while for the U{1) 
gauge symmetry, ki need not be an integer. All the other anomaly coefficients such as 
GiGjGk and [t/(l)^]^ x Gi should vanish. 

In our models, the gauge symmetry, which we are interested in, is SU{3)c x SU{2)l x 
U{1)y X U{1)a- So, the relevant Green-Schwarz anomaly cancellation conditions are 

where ^3 and A2 are the [SU{?>)cf' x U{1)a and \SU{2)l\'^ x U{1)a anomaly coefficients. 
In this paper, we do not consider the anomalies involving the [/(l)y, because its associated 
Kac-Moody level ki is not an integer in general and this condition is not very useful from 
an effective low energy theory point of view 41j. Similarly, the [[/(l)^]^ anomaly can be 
cancelled by the Green-Schwarz mechanism, but this condition also has an arbitrariness 
from the normalization of U{1)a- And the [?7(1)a]^ x U{1)y anomaly does not give any 
useful low energy constraint. 
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A. The Supersymmetric DFSZ Axion Model 



We introduce two SM singlet fields 5* and 5" to break the PQ symmetry, and two SM 
singlet fields X and X to generate the dimension-5 operators for the active neutrino masses 
after they are integrated out. The superpotential is 

W^tot = Wo + W,, (27) 

where 

Wo = ylQ^u]H^ + ytM^'^Hd + yt^Ue^H^ + h-^H,H^ + ys^-^^ , (28) 



= yiLiHuX + y'jnj(f)X + yxSXX + '-(f)^ + Xij^UiUj . (29) 



''PI 

Similar to the discussions in the intermediate-scale U{iy model, we assume that the sum 
of the supersymmetry breaking soft masses for S and 5" is negative, i.e., m| + m|, < 0. 
Then, the 5* and S' fields can acquire intermediate-scale VEVs around 10^*^ GeV, which 
gives us the PQ symmetry breaking scale /„. And the term is given by 

^ = ~ 10^ GeV . (30) 

Mpi 

To forbid the other renormalizable operators and the dangerous non-renormalizable op- 
erators in the superpotential which are allowed by the gauge symmetry and can destabilize 
the axion solutions to the strong CP problem, we introduce a Z102 discrete symmetry arising 
from the breaking of an anomalous U{1)a gauge symmetry. Under the U{1)pq smmetry and 
Z102 discrete symmetry, the charges for the particles in this model are given in Table ITTTl Be- 
cause the anomaly coefficients ^3 and A2 are equal to 60 and 40, respectively, the anomalies 
can be cancelled by Green-Schwarz mechanism if /ca = 3 and k2 = 2, i.e., A^/k^ = A2/k2. 

In addition, using mathematica code, we have shown that up to dimension-5 operators 
in the superpotential (dimension-6 operators in the Lagrangian), the terms in Eqs. (|28p 
and (j29|) are the only operators which are allowed by the gauge symmetry and Z102 discrete 
symmetry. Moreover, this ^102 symmetry forbids to high orders the dangerous terms of the 
following type in the superpotential 
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(31) 



TABLE III: Under the U{1)pq symmetry and Z102 discrete symmetry, the charges for the particles 
in the supersymmetric DFSZ axion model. 





Q^ 


«f 


dt 






m 


Hu 


Hd 


5 


S' 


X 


X 




Uil)pQ 








-1 


-5/4 


1/4 


3/4 





1 


-1/2 


1/2 


5/4 


-3/4 





Z102 


23 


85 


11 


73 


63 


21 


96 


68 


20 


31 


35 


47 


34 



which can potentially destabilize the axion solutions [37[. However, as pointed out in 
Ref. |38i], the axion solutions to the strong CP problem may be destabilized by the non- 
renormalizable terms in the Kahler potential. Therefore, how to stabilize the axion solutions 
from the non-renormalizable terms in the Kahler potential is still an interesting question 
which deserves further study. 

After the PQ symmetry breaking, we obtain the superpotential for neutrino sector 



Wu = XijLiUjHu 



— — + -( 

Mx 3 



(32) 



where 



A 



Viy'j , Mx = yx{S) , m. 



(33) 



PI 



Thus, we obtain the superpotential in the MSSNM. 



B. The Supersymmetric KSVZ Axion Model 

In addition to those particles in above DFSZ axion model, we introduce one pair of 
vector-like fields \1/ and \l/ which belong to the 5 and 5 representations, respectively, in the 
SU{5) language. The superpotential is 

Wtot = Wo + W, + y^S'^^ , (34) 

where Wo and Wi, are given in Eqs. and (f^ . respectively. 

Similar to the DFSZ axion model, to forbid the other renormalizable operators and the 
dangerous non-renormalizable operators in the superpotential which are allowed by the gauge 
symmetry and can destabilize the axion solutions to the strong CP problem, we introduce a 
ZiQ2 discrete symmetry arising from the breaking of an anomalous U{1)a gauge symmetry. 
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TABLE IV: Under the U{1)pq symmetry and Z102 discrete symmetry, the charges for the particles 
in the supersymmetric KSVZ axion model. 





Q^ 




dt 


Li 




m 


Hu 


Hd 


5 


S' 


X 


X 





^ + 


PQ 





























-1 











1 


Z102 


23 


85 


11 


73 


63 


21 


96 


68 


20 


31 


35 


47 


34 


71 



Under the U{1)pq symmetry and Z102 discrete symmetry, the charges for the particles in 
this model are given in Table HVl Note that the non-renormalizable term ys{SS'y/Mpi in 
the superpotential, which violates the U{1)pq symmetry, can be generated from quantum 
gravity interaction. 

Because of the additional fields ^ and "if, we obtain that the anomaly coefficients ^3 and 
A2 are equal to 80 and 60, respectively. And the anomalies can be cancelled by Green- 
Schwarz mechanism by choosing = 4 and k2 = 3. 

The rest discussions are similar to those in above subsection, so, we will not present them 
here. 

VI. CONCLUSIONS AND DISCUSSIONS 

We proposed the minimal supersymmetric sterile neutrino model with late time phase 
transition. The masses for the sterile neutrinos are about 1 eV, and the active neutrino 
masses and the mixings among the active and sterile neutrinos are generated during late 
time phase transition. We can also forbid the dangerous operators by introducing x Z2 
discrete symmetry. In the MSSNM, the current neutrino data from all the experiments 
include the LSND can be explained simultaneously, and one can automatically evade the 
constraints on sterile neutrinos from the BBN, large scale structure surveys and WMAP. 
However, how to produce the 1 eV masses for the sterile neutrinos is an interesting question. 
Also, the supersymmetry breaking is mediated by gauge interactions, then, the fi problem 
is still a severe problem in the MSSNM. To realize the MSSNM naturally, we considered 
the supersymmetric intermediate-scale U{1)' model, the supersymmetric low energy f/(l)' 
model with a secluded f/(l)'-breaking sector, and the supersymmetric DFSZ and KSVZ 
axion models. In these models, the 1 eV masses for the sterile neutrinos can be obtained 
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via the high-dimensional operators by integrating out the heavy fields, and the dimension- 
5 operators for the active neutrino masses and the mixings among the active and sterile 
neutrinos can also be generated by integrating out the heavy fields. Moreover, the /i problem 
can be solved elegantly. Furthermore, for the low energy f/(l)' model with a secluded 
f/(l) '-breaking sector, we briefiy gave a scenario where the sterile neutrino masses are also 
generated during late time phase transition. 
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